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Abstract. It is by now well known that the relativistic heavy-ion collisions at RHIC, BNL have 
produced a strongly interacting fluid with remarkable properties, among them the lowest ever ob- 
served ratio of the coefficient of shear viscosity to entropy density. Arguments based on ideas from 
the String Theory, in particular the AdS/CFT correspondence, led to the conjecture — now known 
to be violated — that there is an absolute lower limit l/47r on the value of this ratio. Causal vis- 
cous hydrodynamics calculations together with the RHIC data have put an upper limit on this ratio, a 
small multiple of 1 /An, in the relevant temperature regime. Less well-determined is the ratio of the 
coefficient of bulk viscosity to entropy density. These transport coefficients have also been studied 
nonperturbatively in the lattice QCD framework, and perturbatively in the limit of high-temperature 
QCI). Another interesting transport coefficient is the coefficient of diffusion which is also being stud- 
ied in this context. I review some of these recent developments and then discuss the opportunities 
presented by the anticipated LHC data, for the general nuclear physics audience. 
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1. Introduction 

Transport coefficients of the QCD matter are of fundamental importance not only because 
they represent an important aspect of QCD, but also because they can be calculated from 
first principles. Trying to extract these coefficients reliably from experimental data and 
evaluating them in various theoretical approaches is a very active area of research today. 

This review is addressed to the general Nuclear Physics audience, a majority of whom 
do not work in the area of relativistic heavy-ion collisions, but have probably heard the 
claim that the RHIC experiments have produced the most perfect fluid ever observed. Be- 
fore I explain the meaning of this claim and present the experimental evidence for it, a few 
introductory remarks are in order. It may be recalled that the phase diagram (pressure vs 
temperature) of the most familiar liquid, namely water is known for long with good accu- 
racy. In particular, the coordinates of the triple point and the critical point are known to 
several significant places, and the various phase co-existence lines are well-determined. In 
contrast, the phase diagram of the strongly-interacting matter or the QCD phase diagram 
(temperature vs the net baryon number density or equivalently the baryon number chemical 
potential) is known only schematically from the experimentalist's point of view. Indeed, it 
is not even known with certainty what are the various phases that occur at high densities. 
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The big idea is to map out quantitatively the QCD phase diagram with the relativistic 
heavy-ion colhsions as an experimental tool. Such experiments have been performed at 
SPS (CERN), are being performed at RHIC (BNL), and will soon be performed at LHC 
(CERN), at successively higher energies: up to -v/s^vAr = 1^, 200 and 5400 GeV, at the 
above three facilities, respectively. It is also necessary to systematically scan the energy 
range up to 200 GeV in order to study the QCD matter at high baryon number density and 
to locate the critical point and the phase transition line predicted by some theories. This is 
being done at RHIC and plans are afoot to do it at FAIR (GSI) and NICA (JINR). 

In nonrelativistic fluid dynamics, the kinematic viscosity {v) is defined as ~ rj/ p 
where 77 is the coefficient of shear viscosity or the dynamic viscosity and p is the density 
of the fluid. It allows us to compare the viscosities of fluids with different densities. (Inter- 
estingly, under standard conditions, water has a lower v than air, although its 77 is higher) 
The dimensionless ratio 77/s serves as the relativistic analog of the kinematic viscosity, s 
being the entropy density.^ (The dimensionless ratio f]ln where n is the number density is 
of no use here because n is ill-defined in relativistic heavy-ion collisions.) 
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Figure 1. Comparison of different fluids with each other. See the text for details. The 
lines are drawn to guide the eye. Figure from Ref. [1]. 

Figure I shows 77/5 as a function of (T — To)/ro where T is the temperature and To the 
critical temperature, for various fluids, namely water, nitrogen, helium and the QCD matter. 
The upper three curves, drawn at the respective critical pressures, exhibit a cusp at Tq. The 
liquid and gaseous phases behave differently because the momentum transport mechanisms 
are different in the two cases; see, e.g., [2]. The points labelled Me5o« Gas have large (~ 50 
%) errors and are obtained from chiral perturbation theory. The points labelled QGP are 
from lattice QCD simulations of [3]. The point labelled RHIC is discussed below. This 



^Strictly speaking, relativistic analog of p is (e + P) = sT where e is the energy density and P the 
pressure. 
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comparison of various fluids explains the statement that the RHIC fluid is the most perfect 
fluid ever observed. 

The point labelled RHIC in Fig. 1 was obtained by matching the eUiptic flow data at 
RHIC with the results of viscous hydrodynamic calculations: Consider a non-central col- 
lision of two identical spherical nuclei as shown in Fig. 2. The nuclei travel parallel to the 
z axis, xy plane is the azimuthal or transverse plane and xz plane is called the reaction 
plane. The overlap zone is shown as the shaded area in the figure. In an ultrarelativistic 
colUsion, the nucleons in the non-overlapping zones continue to travel more or less along 
their pre-collision trajectories, leaving behind the almond-shaped overlap zone. The inter- 
esting observables are governed mostly by the overlap zone which has a very high initial 
energy density. The spatial anisotropy of the overlap zone ensures anisotropic pressure 
gradients in the xy plane. This leads to a final state characterized by momentum anisotropy 
and anisotropic distribution of particles in the xy plane. 
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Figure 2. Non-central collision of two nuclei. 
The triple differential invariant distribution of particles emitted in the final state is given 

by 

^ (fN _ d^N _ d'^N 1 
d^p pxdpxdyd^ pxdpxdy 27r 

where pT is the transverse momentum, y the rapidity and cj) the azimuthal angle of an emit- 
ted particle. The azimuthal distribution is Fourier-decomposed, and the leading coefficients 
vi and V2 are called the directed and elhptic flow, respectively. They provide a measure 
of the anisotropy of the flow in the transverse plane, mentioned above. The importance of 
V2 lies in the fact that it is a measure of pressure at early times and hence the measure of 
thermalization of the quark-gluon matter produced in heavy-ion collisions. 

Hydrodynamics is an effective theory that describes the slow, long- wavelength motion of 
a fluid close to equilibrium. It is a powerful technique because given the initial conditions 
and only the equation of state (EoS) of the matter, it predicts the space-time evolution of 
the fluid. Its limitation is that it is applicable at or near (local) thermodynamic equiUbrium 
only. Hydrodynamics plays a central role in modehng relativistic heavy-ion colUsions. 
State-of-the-art calculations for the RHIC data are based on the relativistic causal dissipa- 
tive hydrodynamics. However, the initial calculations were done in the ideal (t//.s = 0) 
hydrodynamics framework. As an example, see Fig. 3. The broad agreement between 
the data and these initial calculations, in particular the mass ordering of V2{pt), led to the 
claim of formation of an ideal fluid at RHIC. 
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Vsm^i = 200GeV ^^^Au + ^^^Au at RHIC 
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Figure 3. Success of ideal hydrodynamics: Minimum-bias elliptic flow data compared 
with ideal hydrodynamics calculations of Huovinen. Figure from Ref. [4], 

2. Transport Coefficients 

Consider a fluid in equilibrium. If it is perturbed so that the density is no longer uniform, it 
responds by setting up currents which tend to restore the equilibrium. In the linear response 
theory, the current or flux is proportional to the force which is negative of the gradient of 
the density: pu = — DVp. The constant of proportionality (D) is called the coefficient of 
diffusion. Other familiar examples are Ohm's law J = crE = — crV</>, Fourier's law of heat 
conduction Q = — /tVT, etc. A slightly more complicated example involves the transport 
of momentum in response to velocity gradients in an anisotropic medium, the constant 
of proportionality being the shear viscosity tensor Such equations are called constitutive 
equations because they express physical properties of the material concerned. They relate 
the fluxes with the forces, the constants of proportionality being the transport coefficients. 

In addition to the variables such as the hydrodynamic four-velocity, pressure, energy 
density, conserved-number density, etc., hydrodynamics equations also contain transport 
coefficients (shear and bulk viscosities, thermal conductivity, relaxation times, etc.). These 
are external parameters which can be calculated in a variety of ways (see Table 1) and fed 
into the hydrodynamics equations. 

2.1 Transport coefficients from high-temperature QCD 

High-temperature QCD assumes T ^ ^qcd- This is a weak-coupling regime and the 
shear (77) and bulk (Q viscosities can be calculated in the kinetic theory [5], [6], [7]: 
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Table 1. Various ways to calculate transport coefficients for relativistic fluids 



Weak-coupling regime 


Kinetic Theory 


Boltzmann Equation 


Linear-Response Theory 


Kubo I'ormula 


N=4 Supersyimiielric Yang-Mills Theory 




Strong-coupling regime 


Lattice Gauge Theory 


Kubo Formula 


N=4 Supersymmetric Yang-Mills Theory 





It is clear that as the temperature T rises, rj/T^ increases while C,/T^ decreases, and the 
ratio C/r] C(i decreases. Note that the bulk viscosity vanishes for any conformal field 
theory, and QCD becomes conformal in the limit of high T. On the other hand, when 
T ~ 200 MeV, QCD is far from weakly coupled and the above results can provide only a 
rough estimate of r] and (. 

2.2 Transport coefficients from string theory 

Like ordinary fluids black holes too are thermal systems having notions of temperature and 
entropy. An object falling on the surface of a fluid in equilibrium generates disturbance 
which dies down due to dissipative nature of the fluid. Similarly, the black hole horizon 
gets deformed when an object falls on it. However, it soon recovers its equilibrium shape. 
Thus the notion of "viscosity" is applicable to a black hole as well, and the connection 
between hydrodynamics and black-hole physics does not seem very far-fetched. 

Anti-deSitter/conformal field theory (AdS/CFT) correspondence refers to the equiva- 
lence or duaUty between string theory defined on a certain AdS space and a CFT defined 
on its boundary. It allows the calculation of properties of a strongly-coupled CFT in terms 
of those of a weakly-coupled string theory. For an elementary introduction to these ideas, 
see [8]. 

Ordinary quantum mechanics rules out vanishing rj/s for weakly coupled theories, i.e., 
for theories with well-defined quasiparticles [9]. Kovtun et al. [10], using string theory 
methods, showed that r]/s = l/Aw, for a large class of strongly interacting quantum field 
theories whose dual description involves black holes in the Anti-deSitter (AdS) space. The 
value 1 /47r was conjectured to be its absolute lower bound (KSS bound) for aU substances. 
However, it has recently been reaUzed that the KSS bound is violated for certain conformal 
field theories (CFT) [11], [12]. 

Of course, applying above ideas to the fluid produced at RHIC is speculative: QCD near 
the deconfinement transition temperature is not a CFT, and its gravity dual (if it exists) is 
not known. 



2.3 Transport coefficients from lattice QCD 

Figure 4 shows the results of a quenched QCD (no dynamical quarks or niq oc) calcula- 
tion of rj/s from Ref. [3] in comparison with the high-T QCD results quoted in subsection 
2. 1 and the KSS bound mentioned in the previous subsection. However, these initial results 
have now been superseded by more recent results described next. 
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Figure 4. First lattice QCD results onr)/ s. Figure from Ref. [3]. 



Recently Meyer has calculated both 77/5 and s on the lattice, with higher statistical 
accuracy and a more efficient algorithm, assuming SU(3) gluodynamics [13], [14]. He gets 



77/s = 



_ / 0.134(33) (T = 1.65Tc) 



0.102(56) (T=1.24Tc), 



where the errors contain an estimate of the systematic uncertainty. This is consistent with 
the KSS bound. Further, 



0.008(7) 
0.065(17) 



0.15 


" 0.37 
0.01 



(T = 1.65Te) 

(r= 1.24Te), 



C/s = 0.73(3) 



2.0 
0.5 



(T = i.02re), 



where the statistical error is given and the square bracket specifies conservative upper and 
lower bounds. Note the sharp rise in C,/ s just above Tc- Similar dramatic rise was seen 
in the results presented in [15]; see Figure 5. They extracted the bulk viscosity in the 
presence of light quarks by combining low-energy theorems with lattice data on the QCD 
EoS. However, it is now realized that a determination of C, from correlation functions of 
the energy-momentum tensor is more subtle. Shortcomings of this calculation have been 
pointed out in [16]. If the sharp rise of C/s just above Tc is confirmed, it would imply that 
the QGP is not a perfect fluid near Tg! 

For a review of the progress made in extracting transport properties of the gluonic plasma 
from lattice simulations, see [17]. 



2.4 Transport coefficients extracted from RHIC data 

Here the basic idea is that the shear viscosity reduces the elliptic flow: V2 (viscous fluid) < 
'!;2 (ideal fluid). This is easy to understand: recall that V2 is a measure of the flow 
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Figure 5. Bulk viscosity vs temperature. Tiie tiu'ee curves (top to bottom) correspond 
to ujit (T) = 0.5, 1,1.5 GeV respectively, where cjq [T) is a scale at which the perturba- 
tion theory becomes valid. The three curves reflect the uncertainty in the determination 
of this scale parameter. Figure from Ref. [15]. See, however, the discussion in the text. 

anisotropy in the azimuthal plane. Viscosity is the result of a frictional force. Frictional 
force being proportional to the flow velocity has a relatively stronger effect on fast-moving 
particles emerging in the reaction plane. This reduces the anisotropy and hence V2. 

Thus if one has a good control on V2 (ideal fluid), one can adjust r//s to fit the data on V2, 
and thus extract 77/5. This has been done by several groups [18], [19], [20], [21], [22], [23], 
[24]. Results of [20] are shown in Fig. 6. The Glauber and colour-glass-condensate (CGC) 
initial conditions for hydrodynamic evolution yield 77/s ~ 0.08 and ~ 0.16 respectively, 
showing the sensitivity of the extracted rj/s to the initial conditions. 



Glauber CGC 




Figure 6. Minimum-bias elliptic flow data for charged hadrons in 200 GeV Au-Au 
collisions; only statistical errors are shown. Line thickness is an estimate of the accu- 
mulated numerical error. Figure from Ref. [20]. 



3. Diffusion 

Transport of heat occurs by diffusion (Fourier's law) as well as by propagating waves 
(Maxwell-Cattaneo law). The latter mechanism underlies the well-known phenomenon 
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of second sound in superfluid helium (Table 2). Similarly the transport of momentum 
too occurs by diffusion (Navier-Stokes equation) or by propagating waves. These are the 
transverse or shear waves^, different from the more familiar longitudinal or compressional 
waves. For a detailed discussion see [25]. 

Table 2. Various transport phenomena 



Transport of 


By diffusion 


By propagating waves 


Experimental situation 


Heat 


Fourier's law 


Maxwell-Cattaneo law 


Second sound in superfluid He 


Momentum 


Navier-Stokes eq. 


Maxwell-Cattaneo law 


Propagating shear waves 


Conserved no. 
e.g.,B,Q,S 


Pick's law 


Kelly's law 


RHIC? 



We recently studied the transport of a conserved number such as the net baryon number 
{B), charge (Q), or strangeness (S) conserved in strong interactions, in the acausal and 
causal, or the first- and second-order theories of relativistic diffusion [26]. We found that 
Pick's diffusion smooths out gradients in the number density monotonically. In contrast, 
in Kelly's theory the gradients may be transiently amplified, i.e., the density profile may 
stiffen at intermediate times; see Fig. 7. We proposed experimental observables and argued 
that the RHIC data can potentially distinguish between the above two mechanisms. If the 
second-order theory of diffusion is ruled out by the data, one gets a handle on the relaxation 
time and hence on thermalization. 




Figure 7. Conserved number density profile as a function of time and space-time 
rapidity in the Pick (left) and Kelly (right) theories of diffusion. Pigure from Ref. [26]. 



^They are used, e.g., in Magnetic Resonance Elastography (MRE) for quantitatively imaging ma- 
terial properties. 
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4. Viscosity in Ordinary Finite Nuclei 

This review is about the transport properties of the fluid produced at RHIC. But it is perti- 
nent to ask what do we know about the transport coefficients of the matter in an ordinary 
finite nucleus at low temperatures. Hydrodynamic models have a long history in Nuclear 
Physics — recall the success of the Uquid drop model. The coefficient of shear viscosity of 
a finite nucleus can be obtained from (i) analysis of the widths of giant resonances within 
the hydrodynamic model, (ii) the process of fission studied within the liquid drop model, 
and (iii) kinetic theory. This has been discussed recently in [27] . Using entropy density for 
a free Fermi gas or for noninteracting nucleons in a Woods-Saxon potential, they obtained 
values of ri/s larger than but not drastically different from those for the RHIC fluid. 

5. What about LHC? 

It is clear from Table 3 that at LHC we would need r]/s and (/s up to T AT^. Some 
preUminary lattice results are now available, but for the SU(3) pure gauge theory (i.e., 
quarkless QCD) [17]. For a careful extraction of r//.s and (/s from LHC (and RHIC) data, 
we need to incorporate the T-dependence of these transport coefficients in hydrodynamic 
calculations, among other refinements of these calculations; see the next section. 

Table 3. Comparison of central Au- Au collisions at RHIC and central Pb-Pb collisions 
at LHC 





RHIC (Au-Au) 


LHC (Pb-Pb) 


V^NN 

Initial temperature 
Initial energy density 
Lifetime 


200 GeV 
~2Tc 

~ 5 GeV/fm^ 
~ 10 fm/c 


5.5 TeV 

~4Tc 

15-60 GeV/fm^ 
> 10 fm/c 



6. Take-Home Message 

• Elliptic flow at RHIC has put a robust upper limit on the value of 77/5 of the RHIC 
fluid: r]/s < 5/ (47r). This is the average value in the relevant temperature region. 

• Uncertainties associated with the initial conditions, EoS, bulk viscosity, hadronic 
stage, freezeout procedure, different versions of the second-order hydrodynamic 
equations prevent a more precise determination of rj/s. See Ref. [28] for a detailed 
discussion. 

• The bulk viscosity is not yet weU-determined in the deconfinement transition region. 

• Analysis of the RHIC data, which might throw Ught on the coefficient of diffusion, 
is awaited. 
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